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SWIARY 

This  paper  presents  the  results  of  a  series 
of  experiments  on  the  traveling— salesman  problem. 
The  purpose  of  these  experiments  was  to  Investigate 
the  efficiency  of  the  lineaivprogrammlng  technique 
In  solving  this  type  of  problem. 
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SOMB  BXPBRIHBNT3  ON  THB  TRAVELING— SAIJS3MAH  PROBLEM 

by 

J.  T.  Robacker 

1.  INTRODUCTION  AND  DISCUSSION  OF  RESULTS 

During  th«  early  part  of  195* ,  Dantzig,  Fulkerson,  and 
Johnson  formulated  a  linear-programming  teohnique  for  the 
solution  of  the  traveling-salesman  problem  [l]  .  The  methods 
employed  by  these  authors  quickly  generated  optimal  tours 
when  a  good  tour  was  known.  The  question  naturally  arose 
concerning  the  efficiency  of  the  technique  when  a  random 
tour  was  chosen  as  the  starting  point.  Dr.  I.  Heller  con¬ 
jectured  that  it  might  be  efficient  and  proposed  that  experi¬ 
ments  be  made  as  one  means  to  investigate  this  question  and 
to  gain  a  greater  degree  of  insight  into  the  dynamics  of  the 
technique  itself. 

A  series  of  10  nine-city  traveling-salesman  problems 
was  considered.  In  each  case  the  distances  between  cities 
were  chosen  from  a  table  of  two-digit  random  numbers.  The 
starting  tour  was  then  chosen  to  be  in  the  order  of  the 
natural  numbers  (with  the  cities  numbered  1  through  9  at 
the  outset).  In  the  last  problem,  the  technique  was  applied 
first  to  find  the  maximal  tour  and  then  to  proceed  from  it 
to  the  minimal  tour. 
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Tables  1  to  10  are  the  distance  tables  for  the  10 
problems.  An  element  ajj(aij  “  aji^  ln  distance  table 
represents  the  distance  from  the  i— th  city  to  the  th  city. 
It  is  to  be  noted  that  since  these  distances  were  chosen 
from  a  table  of  random  numbers,  they  cannot  be  construed 
necessarily  to  represent  straight— line  distances  in  Euclidean 
space.  At  the  foot  of  each  distance  table  is  a  table  listing 
the  tours  obtained  at  the  end  of  each  iteration*  (the  last 
entry  being  the  optimal  tour).  Table  11  lists  the  tours 
obtained  by  starting  with  the  maximal  tour  associated  with 
Table  10  and  proceeding  to  the  minimal  tour. 

The  largest  number  of  iterations  needed  was  six,  while 
the  average  was  only  a  little  under  four.  The  average  Mme 
to  work  one  of  the  problems  was  about  5  hours.  Since  for 
a  nine— city  problem  there  are  ~  ■  20,160  possible  tours, 
it  is  apparent  that  the  simplex  method  which  was  employed 
was  extremely  efficient.  In  addition,  it  is  noteworthy  that 
the  only  secondary  constraints  introduced  in  the  solution  of 
these  problems  were  upper  bounds. 

In  connection  with  these  experiments,  A.  W.  Boldyreff 
suggested  an  approximation  procedure,  the  merit  of  which 

♦One  iteration  is  defined  to  be  the  process  of  going  from 
one  tour  to  another  tour  of  the  same  or  shorter  length. 
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11*8  In  Its  Inherent  simplicity  and  In  the  rapidity  with 
which  it  may  be  applied.  An  application  of  this  approximation 
method  to  the  49-city  problem  of  [Y]  gave  a  tour  of  851  units 
as  compared  with  the  optimal  of  699  units,  an  error  of  21$6. 

2.  BRIEF  REVIEW  OF  THE  TECHNIQUE  OP  SOLUTION 

The  mathematical  technique  found  in  [Y]  will  briefly 
be  reviewed  here,  and  will  then  be  illustrated  by  way  of  an 
example . 

The  factors  ,  referred  to  as  "potentials,"  are  com¬ 
puted  from  the  formula 

(l)  Ti  +  ~  alJ  *  0  for  a11  (i»j)  ln  the  basis, 

i.e.,  for  all  (i,j)  corresponding  to  column  vectors  in  the 
basis. 

Computationally,  this  means  first  finding  a  loop  with 
an  odd  number  of  links  which  correspond  to  vectors  in  the 
basis;  then  starting  at  a  city  k  of  the  loop,  adding  and  sub¬ 
tracting  the  lengths  of  the  links  alternately;  the  resulting 
number  is  then  2  t^.  Having  Tk,  we  can  then  compute  all  the 
from  ( l) . 
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Figure  1 

We  next  choose  the  link  (k,/)  such  that 
(2)  Tk  +  ri  ~  V<  -  n»*  (»*  +  Tj  -  aij>  >  0 

and  allow  the  value  of  the  corresponding  variable  xkj  to  be  equal to  0 . 
The  values  of  the  basic  variables  must  now  be  corrected  so 
that  the  sum  of  the  values  of  link  variables  around  any  city 
must  equal  2.  For  this  purpose  we  try  to  find  an  even  loop 
containing  (k ,1).  For  example,  suppose  the  starting  tour 
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la  illustrated  aa  i.i  Figure  1  along  with  the  potentials  and 
links  corresponding  to  basic  variables.  Hie  numbers  on  the 
links  are  the  lengths  of  those  links  as  chosen  from  Table  1. 

It  is  seen  that  the  greatest  value  which  (6,2)  may  have  is 
0  since  the  value  of  each  link  in  the  basis  is  bounded  above 
by  1.  Consequently,  we  put  a  bar  on  (9,1 ),  indicating  that 
it  has  attained  its  upper  bound  and  is  replaced  in  the  basis 
by  (6,2).  Since  (9,1 )  is  no  longer  in  the  basis,  we  must 
recompute  the  prices  according  to  (l). 

Figures  1—5  illustrate  the  stepB  necessary  to  complete 
the  first  iteration.  The  solid  lines  are  those  in  the  tour. 

In  Pigure  5  a  new  tour  is  achieved.  By  repeating  the  process 
we  ultimately  reach  the  optimal  tour  when 

n  +  Tj  -  a«  * 0 


(3) 


for  all  i, J  -  1,2, . . . ,9. 


© 
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Flgure  5 


-15  i 


! 


Figure  4 
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AH  APPROXIMATION  TECHNIQUE 

The  following  approximation  procedure  has  been  suggested 
by  A.  W.  Boldyreff. 

Start  with  the  first  three  cities.  There  is  then  on3y 
one  tour  of  these,  namely,  125.  We  no"  8efc 


P(l,J;^)  •  ai4  ♦  aj4  ”  aiJ 


i,J  -  1,2,5 


Let  k  and  i  be  the  cities  such  that 


(5)  F(k,J;4)  -  rain  F(i,J;,0  • 

i»J"l*2*5 

Then  our  tour  of  the  first  four  cities  replaces  the  link 
(k,|)  in  the  tour  of  the  first  three  cities  by  the  two  links 
(k, 4)  and  (4,1).  We  repeat  this  process,  adding  one  city 
at  a  time,  until  we  have  a  tour  of  all  the  cities. 

To  illustrate  let  us  consider  Table  1  and  tabulate  the 


procedure  as  follows. 


m- 1521 
7-88-55 
-9- 


12  3  1  city 

4  46  98  length  of  link 


longth  of  link 
*15 


26 

37 

46 

26 

*14 

1 

2  * 

3 

4 

1 

4 

46 

46 

26 

84 

50 

95 

54 

84 

1  2 

5 

*  3 

4 

1 

4 

50 

95 

46 

26 

70  19 

79 

46 

96 

70 

.  X 

1  2 

5 

6 

3 

4 

4  50 

79 

46 

46 

26 

longth  of  link 


77 

88 

66 

75 

26 

19 

*17 

5  X 

6 

3 

4  7 

1  city 

50 

79 

46 

46 

26 

19 

longth 

74 

98 

16 

62  36 

76 

55 

74 

_  X 

1 

2 

5 

8  6 

3 

4 

7 

4 

50 

16 

62  46 

46 

26 

19 

longth  of  link 


o  c  -*  k  7  9  1  tour 

1  4  2  50  16  62  46  46  26  29  2  length  of  ling 


At  each  etage  the  X  mark.  the  link  which  la  replaced. 

In  thl.  example  the  final  tour  ha.  a  length  of  281  while  the 
optimal  tour  ha.  a  length  of  252.  The  error  la  21*.  In  moat 
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cases  considered  here  the  error 


is  less  than  this,  as  can  be 


seen  from  the  comparison 
There  are  many  ways 


in  the  tables  of  Section  5- 
In  which  the  accuracy  of  this  tech- 


nique  may  be  Improved,  but  it  is  felt  that  the  vaulting 
complexity  would  probably  re.ult  in  a  relatively  email  increaae 


in  accuracy. 


City 
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i 


Iteration 

Tour 

Length  of  Tour 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

432 

1 

1 

2 

6 

3 

4 

5 

8 

7 

9 

271 

2 

1 

2 

6 

3 

8 

5 

4 

7 

9 

232 

Approximation 

Technique 

1 

2 

5 

8 

6 

3 

4 

7 

9 

281 

Table  1 


City 
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Iteration 

Tour 

Length  of  Tour 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

459 

1 

1 

2 

3 

4 

5 

6 

7 

9 

8 

382 

2 

1 

3 

4 

5 

2 

6 

7 

9 

8 

336 

3 

1 

3 

4 

7 

6 

5 

2 

9 

8 

240 

4 

1 

3 

4 

7 

6 

8 

9 

2 

5 

232 

5 

1 

3 

4 

5 

2 

9 

8 

6 

7 

220 

6 

1 

3 

4 

7 

5 

2 

9 

8 

6 

204 

Approximation 

Technlaue 

1 

4 

7 

5 

2 

9 

8 

6 

1 

225 

Table  2 


95 


63  17 


23  69  60 


21  72 


55  49  20  86  75 

79  32  50  57  23  53 

42  80  99  61  29  16  78 

98  73  48  8  15  2  33 


Iteration 


Tour 


1  4  8  6  3  2  5  7  9 


14  7  952368 


Length  of  Tour 


L  2 

3 

4 

5 

6 

7 

8 

9 

L  4 

5 

2 

3 

7 

6 

8 

9i 

1  4 

3 

2 

5 

7 

6 

8 

9 
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2  63 


3  64  95 


4  24  18  76 


£5  76  85  74  15 

_ _ 

6  38  65  68  76  11 


7  38  90  90  75  13  31 

8  51  43  77  92  26  6 

9  92  47  54  55  8l  32 

1  2  3  4  5  6 

City 


Iteration  Tour 

0  123456789 

~1  1298765  .4 _ 3_ 

2  1245678 _ 2 — 1_ 

3  124  57689  3- 

”4  139245867L 

5  1  3  9  6  8  2  4  5  7. 

Approximation  186754293 
Technique _ _ — - — - 

Table  5 


Length  of  Tour 


Length  of  Tour 


RM-1521 


BS 


KM— 1521 
7-28-55 
-22- 


R3FKRBKCE3 


W 


Dantzlg,  0.,  D.  Fulkerson,  and  S.  Johnson, 
of  a  Largo  Scalo  Travollng  Salesman  Problem 
of  the  Ooer.  Res.  Soc.  of  Araer.,  Vol.  2,  No 
Nov.,  1554,  pp" 


Solution 
"  Jour. 

*7 


UNCLASS/F/E 


NMtoSlWVSlIWUll 


Reproduced  by 

DOCUMENT  SERVICE  CENTER 

KMTT  BUILDING,  DAYTON,  2,  OHIO 


This  document  to  (J 
ration  of  the  contract  a. 
to  the  folknrlnf  Mhdr 
Docooeot  Service 


property  of  the  United  States  Government.  It  is  furnished  for  the  du- 
I  shall  be  returned  when  no  longer  required,  or  upon  recall  by  ASTIA 
>•:  Armed  Services  Technical  Information  Agency, 
Center,  Knott  Building,  Dayton  2,  Ohio, 


NOTICE  WHEN  QOS  IftKMXNT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
Anr  USED  for  any  PURPOSE  OTHER  than  in  connection  with  a  DEFINITELY  RELATED 
COVERS  VOL  NT  PROCUPUEVCKNT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPOKSimurY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  RAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
RAID  DRAWINGS  f FKCtnCATIONS.  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  0T« I  ftWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER. 
PERSON  OR  COnrOBUTXJfv  OH  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  M  ANUFACTURE 
USE  OR  SELL  ANY  PJ.T  t  NT  ED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


